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Abstract 

Abstract The study of the properties of cosmic structures in the universe is one 
of the most fascinating subject of the modern cosmology research. Far from being 
predicted, the large scale structure of the matter distribution is a very recent discov- 
ery, which continuosly exhibits new features and issues. We have faced such topic 
along two directions; from one side we have studied the correlation properties of the 
cosmic structures, that we have found substantially different from the commonly 
accepted ones. ^From the other side, we have studied the statistical properties of 
the very simplified system, in the attempt to capture the essential ingredients of 
the formation of the observed strucures. 



1. INTRODUCTION 

The existence of a cosmic structures in the universe is one of the most fascinating findings 
of the two last decades in observational cosmology. Galaxy distribution is far from homo- 
geneous on small scale and large scale structures (filaments and walls) appear to be limited 
only by the boundary of the sample in which they are detected. There is currently an acute 
debate on the result of the statistical analysis of large scale features. In the past years 
we propose a new statistical approach, which has shown surprisingly a fractal structure, 
extending from small scales to 

distance beyond to AOh~^ Mpc and 
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even up to lOOh Mpc. In the following wc report our findings along two directions of 
investigations: in the first part we describe the standard analysis and its limits of applica- 
bility. We describe our novel analysis and the results we got in the characterization of large 
scale structures. Contrary to standard claims, this analysis is completely consistent with 
the results of the standard analysis: it's the interpreation of the latter which is radically 
different. In the second part, we refer on a study of the dynamics of a very simple model 
of gravitational formation of structures. We have analysed the evolution and the statistical 
spatial properties of a N-body system of point masses, interacting through gravity. The 
system is arranged as to simulate an infinite system of particles and with very simple initial 
conditions. 



2. COSMIC STRUCTURES 

The usual way to investigate the properties of the galaxy spatial clustering is to measure 
the two point autocorrelation function ^(r) This is the most used statistical tool, since it 
can be measured quite accurately with current redshift surveys, ^(r) is the spatial average of 
the fluctuations in the galaxy number density at distance r, with respect to an homogeneous 
distribution with the same number of galaxies. 

Consider a little volume 6V at position ff, let n(r^) be the density of galaxies in SV and 
< n >= N/V the density of galaxy in the whole sample. 

The galaxy density fluctuations in 6V with respect to the average galaxy density < n >, 
i.e. the galaxy relative density fluctuations, is: 

n{fl)- <n> _ 6n{ri) 
<n> ~ <n> 

The two point correlation function ^(r) at the scale r is the spatial average of the product 
of the relative density fluctuations in two little volumes at distance r: 

'^"'(^"^ + ^) ^'^(^) y _ < n{fi)n{fi + r)>i ^ 
<n><n>*~ <n>2 

where the average is performed over the sample. Roughly speaking, a set of points is 
correlated on scale r if ^(r) > 0; it is uncorrelated over a certain scale R \i ^{r) =0 for 
r > R. In the latter case the points are evenly distributed at scale R > r ot, in another 
words, they have an homogeneous distribution at scale R > r. In the definition of ^(r), 
the use of the sample density < n > as reference value for the fluctuations of galaxies is 
the conceptual assumption that the galaxy distribution is homogeneous at the scale of the 
sample. 

Clearly such an approach is valid if the average density < n > of the sample is the 
average density of the distribution, or, in other words, if the distribution is homogeneous 
on the scale of the sample. For this reason, ^ (r) analysis assumes the homogeneity and it is 
unreliable for testing it. 

In order to use ^(r) analysis, the density of galaxies in the sample must be a good 
estimation of the density of the whole distribution of galaxies. This may either be true or 
not; in any case, it should be checked before applying ^(r) analysis ^. 

In addition to such criticisms, the usual interpretation of ^(r) measure is uncorrect for 
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spatial distribution of points with respect the amplitude of the ^(r). The correlation length 
of the distribution rg is indeed defined as the scale such that '^(rg) = 1 

Such a definition is uncorrect since, in statistical mechanics, the correlation length of the 
distribution is defined by how fast the correlations vanish as a function of the scale, i.e. by 
the functional form of the ^ (r) and not by its amplitude . 

The quantity ro, then, does not concern the scale of fluctuations and it is not correct 
referring to this as a measure of the characteristic size of correlations 10,11^ According to 
the ^(r) definition, ro simply separates a regime of large fluctuations 6n/ < n >S> 1 from 
a regime of small fluctuations 5n/ < n ><C 1: this is correct if the average density of the 
sample < n > is the average density of galaxy distribution . 

Such problems of ^(r) approach can be avoided analysing the spatial correlations of the 
data set without any a priori assumptions on the homogeneity scale of the data itself. ^ 

The way to perform such an unbiased analysis is to study the behaviour of the conditional 
average number of galaxies < N{< r) > or the average conditional galaxy density r*(r) 
versus the scale r. The two quantities are respectively: 



< iV(< r)>=B- r° (3) 

where iV(< r) is the number of galaxies contained in a sphere of radius r centered on 
a galaxy of the sample. < N{r) > is the average of N(< r) computed in all the spheres 
centered on every galaxy of the sample. 

r*(r) is then the corresponding average density of galaxies in spheres of radius r The 
exponent D is called the fractal dimension and characterises in a quantitative way how the 
system fills the space, while the prefactor B depends on the lower cut-off of the distribution. 

< A^(< r) > and r*(r) are the suitable statistical tools to detect the two-point correla- 
tion properties of a spatial distribution of objects and the possible crossover scale between 
different statistical distributions. 

If the point distribution has a crossover to an homogeneity distribution at scale R, T*{r) 
shows a fiattening toward a constant value at such scale. In this case, the fractal dimension 
of Eq.(3) and Eq.(4) has the same value of the dimension of embedding space d, D = d (in 
three-dimensional space D = 3)^^'^'^^. 

If this does not happen, the density sample will not correspond to the density of the 
distribution and it will show correlations up to the sample size. The simplest distribution 
with such properties is a fractal structure A fractal consists of a system in which more 
and more structures appear at smaller and smaller scales and the structures at small scales 
are similar to the ones at large scales. The distribution is then self-similar. It has a value 
of D smaller than d, D < d. In 3-dimensional space, d = 3, a fractal has D < 3 and 

r*(r) is a power law. The value of N{< r) largely fluctuates by both changing the starting 
point, from which we compute A^(< r), and the scale r. This is due to the scale invariant 
feature of a fractal structure, which does not have any characteristic length jg 
simple to show that if we analyse a fractal structure with ^(r), we can obtain a value for 
the correlation length ro, which evidently does not have any relation with the correlation 
properties of the system. In fact such a value is simply a fraction to the size of the sample 
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Figure 1: r*(r) in the range of scales 0.5 lOO/i ^Mpc for all the avalaible 3D galaxy data. 
A reference line with a slope —1 is also shown (i.e. fractal dimension D = 2). 



According to our criticism to the standard analysis, we have performed the measure of 
galaxy conditional average density r*(r) in all the three dimensional catalogs available. Our 
analysis is carried out on several 3D galaxy samples. The results are are collected in Fig. 

1 

r*(r), measured in different catalogues, is a power law as a function of the scale r, 
extending from ~ 0.5 to 30 — AOh~^ Mpc, without any tendency towards homogenization 
(flattening) 12. In a single case, the LEDA sample it is possible to reach larger scales. 
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smaller scales. This data sample has been largely criticised, but to our knowledge, never in 
a quantitative way. The statistical tests we performed show clearly that up to 50h~^ Mpc 
the results are consistent with all other data 

Such results imply that the ^(r) analysis is inappropriate as it describes correlations 
as deviations from an assumed underlaying homogeneity. The galaxy distribution shows 
instead fractal properties at least in the range r 0.5 — 30 — AOh^^Mpc, which seem to 
extend in a single sample up r ~ lOOh^^ Mpc 

By consequence, the correlation length ro, i.e. the amplitude of the ^(r), for samples 
with such a linear extention, should be a fraction of the sample size: ro should be larger for 
samples whose size is larger. 

This is evident in fig. (2), where we plot the results of the standard analysis ^(r) per- 
formed on the same data sets analysed with r*(r) (fig. (1)). 

In this case, then, ro has no relation with the correlation properties of the system; its 
variation in different samples is not related to any variation of the clustering of the corre- 
sponding data set. 

3. SIMULATIONS OF GRAVITATIONAL CLUSTERING 

The study of the formation of the cosmic strucures we have analysed in the above section is 
one of the most challenging problems in astrophysics. Gravity is the most natural candidate 
for the explanation of the variety of structures we observe. Indeed, the range of scales on 
which the gravitational clustering takes place is really impressive: from 1Q~^ pc to W^pc 
{Ipc = 3.2615 light — yr = 3.0856 a; 10^^ cm). This implies interactions of gravity with 
other physical processes depending on the scale: from turbulence in cold molecular clouds 
to cosmological expansion above galaxy cluster scale. Because such a richness of physical 
processes can be involved in modelling the various structures we observe in cosmos, it is 
actually very difficult to retrieve a clear picture of the statistical properties of self-gravitating 
system. Current astrophysical simulations have reached a high level of refinement, both in 
resolution and in the number of different physical processes which they take into account. 
Such characteristics allow them to study in great detail the single physical problem for which 
they are developed ^. On the other hand they don't allow a clarification the common role 
and the peculiarities of gravitational interaction. On the contrary, we have tried to focus 
on such features analysing the most simplest case of a many-body infinite self-gravitating 
system, without any other ingredient but the gravity. The theoretical approach to such a 
system goes back to Newton himself ^, although it has faced by very few authors (e.g. ^). 
Indeed, the current theoretical effort is quite different since it is devoted to the study of 
evolution of a continuos gravitating fluid, which is assumed to have peculiar initial density 
fluctuations ^Prom the point of view of statistical mechanics, it is very hard to study 
the properties of an inflnite system of self-gravitating particles. This is mainly due to the 
long range nature of gravitational potential, which is not shielded by the balance of far away 
charges, as e.g. in a plasma. Therefore all scales contribute to the potential energy of a 
particle. The peculiar form of the gravitational potential produces two classes of problems: 
those due to the short range (i.e. r — > 0) divergence and those due to the long range (i.e. 
r — ^ -|-oo) behaviour. The former is not uncommon, since it is the same problem which arises 
in electromagnetism. The divergence would cause, e.g., the Boltzmann factor to diverge in 
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Figure 2: ^[r) measure in various VL galaxy samples. The general trend is an increase 
of the ^(r) amplitude for brighter and deeper samples. In the insert panel we show the 
dependence of correlation length ro on sample size Rg for all samples. The linear behaviour 
is a consequence of the fractal nature of galaxy distribution in these samples. 
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The physical nature of this cut-off may be due to many effects, e.g. the dynamical emergence 
of angular momentum barriers. The long range behaviour is of much more concern and is, 
in fact, the problem. It is an easy exercise to verify that the energy of a particle in an 
infinite self-gravitating system diverges. This causes the energy to be non- extensive. As a 
consequence, a thermodynamical limit is not achieved, since as the number of particles goes 
to infinity, even keeping the density constant, the energy per particle diverges. Strangely 
enough, such a problem has not been fully appreciated by many physicists in the field (see 
e.g. ^), as they try to avoid the long range divergence by putting the system in a box "as 
it is usually done with ordinary gas" . In fact, the difference is that in ordinary gas, when 
confining the system in a box, the energy per particle is equal to a constant plus a surface 
term that goes to zero in the thermodynamical limit. In self-gravitating systems, due to 
non extensivity, the energy per particle is neither a constant, nor the surface term goes to 
zero (in fact, it is of the same order of magnitude as the potential energy due to particles 
belonging to the system). 

Another very interesting consequence, which is often not appreciated, is that the ther- 
modynamical definition of temperature, as the parameter which controls the equilibrium of 
the system, doesn't hold for a self gravitating system, since one cannot divide a system into 
smaller subsystems with the same thermodynamic properties of the larger system. 

As a consequence of such difficulties, a satisfying thermodinamic equilibrium treatment 
of such systems is still lacking. 

However we are much more interested in what happens out of equilibrium, during the 
evolution of a system. 

The system we intend to simulate is a infinite many body system. At this aim, the 
particles we effectively consider are confined in a cube of size L submitted to periodic 

boundary conditions. Every particle in the simulation box interacts with all other particles 

and with the periodic replicas of the whole system. 
The initial conditions we consider are: 

(1) random (white noise) initial positions of particles; 

(2) no cosmological expansion; 

(3) zero initial velocities; 

(4) equal mass particles; 

Some snapshots of the temporal evolution of the system, with = 32000 particles, are 
shown in Fig. (3). The time evolution goes from the top to the bottom and the initial 
unclustered distribution of mass points evolves toward a clustered distribution. For each 
snapshot, we plot on the right, the corresponding r*(r). A typical time for the evolution of 
the system is r = 1/ -v/Gp, where G is the gravitational constant and p is the density of the 
system. It is roughly the time needed to a particle to cross the system. 

Fig. 3 shows some interesting features, that we summarize in the following table. 

4. CONCLUSIONS 

The standard analysis of the correlation properties of the galaxy distribution is performed 
through the measure of the ^(r) function. The latter can provide the correct information 
z/the set under analysis is homogeneous inside the sample size. For this reason ^(r) is not 

1 • 1 1 p II- 1 • 1 




igure 3: on the left: snapshots of system evolution. The corresponding time is: t - 
= 2/3t, t = t and t = 4 • r on the right the corrisponding T*{r) for the snapshots on 



Scaling in cosmic structures 9 



t time 


Description 


r*(r) 


to = 


the system is composed by Nparticles 
at rest with spatial possonian distribu- 
tion in the simulation box 


the system has a costant number den- 
sity T*{r,to) = 1 at all the scales. 
At small scale, r*(r) is more fluctuat- 
ing, because of the larger poisson noise 
at small scale 


ti 2/3r 


the system starts to cluster at small 
scale 


r*(r, ti) develops a larger amplitude at 
small scale 


f2 ~ T 


the clustering proccs evolves with the 
merging of the small clusters in bigger 
ones 


the shape of T*{r,t2) appears to be 
quite independent from time, but shifts 
toward larger scales for increasing time 




all the clusters have merged in a sin- 
gle big one, with dimension comparable 
with the simulation box size 


r*{r,tf) does not evolve anymore and 
the system has reached a stationary 
state 



Table 1: on the left: projection of the simulation box onto x — y plane at different time, on 
the right: corresponding measure of r*{r,t) at time t. 



This should be checked before the use of (,{r) and it is possible through T*{r) analysis of 
the sample set. Such an analysis has been performed for the available 3D galaxy samples, 
with the result that the galaxy distribution appear fractal from 0.5 to 100h~^ Mpc. 

For such a range of scales, the ^(r) analysis does not give the correct informations on the 
statistical properties of the galaxy distribution. 

To investigate the formation of such a fractal structure in the universe, we are performing 
simulations of an N-body infinite self-gravitating system. We started from the simplest 
initial conditions and analysed the system with the aforementioned T*{r) function. 

Simulations with a different number of particles have shown that the shape of r*{r,t) 
ioT t < T is independent from N. On the contrary, this is not true for the final state of 
equilibrium, when the system has formed a single cluster 

These measures seems to show that the transient phase, 

during which the collapse occurs, posses a well defined 

thermodynamical limit, which we are currently analysing 

The discrete nature of the N-body system seems to be a 

fundamental ingredient in the development 

of the spatial correlations. 

The latter, indeed, grow at the small scale, where the discretness 
of the point distribution has to be taken into account. 
At the moment is not clear if such a system can develop fractal 
correlations as seen in the galaxy distribution. 
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